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Abstract. Microquasars are accreting X-ray binary systems with non-thermal radio jets. In some of these systems 
the jet is expected to be strongly misaligned with the perpendicular to the orbital plane. If the donor star is 
an early-type star, the jet could collide with the stellar wind producing a standing shock between the compact 
object and the stellar surface. Relativistic particles injected by the jet can be re-accelerated and isotropized at 
the colliding region. If the jet has hadronic content, TeV protons will diffuse into the inner, dense wind leading to 
gamma-ray and neutrino production from interactions with the matter of the wind. In the case of very powerful 
jets, the wind pressure can be overbalanced and the jet might impact directly onto the stellar surface. We present 
estimates of the gamma-ray and neutrino luminosities for different sets of parameters in these scenarios and we 
briefly discuss the effects of this radiation on the donor star and its detectability with current instruments. 
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1. Introduction 

Microquasars (MQs) are accreting X-ray binary systems 
that have shown persistent or episodic non-thermal jet-like 
ejections (Mirabel & Rodriguez 1999). It is believed that 
a significant fraction of all accreting black hole (BH) and 
neutron star (NS) X-ray binaries are likely to exhibit MQ 
behavior at some epoch during their lifetime. Although it 
is not observationally required, the jet axis and the or- 
bital plane of a given MQ system are usually assumed 
to be approximately perpendicular. However, Maccarone 
(2002) has shown that the timescale for the BH spin to 
align with the orbital angular momentum in binary sys- 
tems is often longer than the lifetime of the system it- 
self. Additionally, tidal effects upon the disk can induce 
a precession which could result in a severe misalignment 
(e.g. Larwood 1998, Kaufman et al. 2002). Precession and 
strong inclination of the jet can also be the result of the 
radiative warping of the inner accretion disk (e.g. Wijers 
& Pringle 1999, Ogilvie & Dubus 2001) or might be due 
to spin-spin interactions (e.g. Bardeen & Petterson 1975, 
Sarazin et al. 1980). Observationally, there exists sugges- 
tive evidence for jet misalignments with the normal to the 
orbital plane in wind-fed accreting BHs. For instance, the 
jet axis lies at no more than ~ 35 degrees from the orbital 



plane in V4641 Sgr (see Butt et al. 2003 and references 
therein). Other systems like SS 433, LSI +61 303, and 
Cygnus X-3 are known to have precessing jets (see Katz 
1980, Massi et al. 2004, and Mioduszewski et al. 2001, 
respectively). Precessing jets have been also suggested in 
the case of Cygnus X-l (Romero et al. 2002). 

In this paper we will consider the situation of a mis- 
aligned MQs with its jet lying close to the binary plane 
and being periodically directed straight to the companion 
star. Such a situation seems not to be unlikely. Butt et 
al. (2003), for instance, have considered the induced nu- 
cleosynthesis that should occur on the stellar surface as a 
consequence of the jet impact in such systems. Here, we 
will show that in the case of an early-type companion star, 
the jet, under some conditions, could not reach the stellar 
surface since it is balanced by the strong stellar wind. 
However, another interesting feature appears: a strong 
shock is formed between the compact object and the star 
and protons injected by the jet can be re-accelerated and 
isotropized there. The most energetic of these protons then 
diffuse through the inner wind producing gamma-rays and 
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neutrinos through p 
reaction channels. 

On the other hand, very powerful jets can overcome the 
pressure of the wind and directly impact on the star. In 
such a case the jet energy is dissipated through hadronic 
and electromagnetic cascades that develop in the stellar 
atmosphere. No continuum gamma-ray source above MeV 
energies is expected in this case, but a significant flux of 
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high-energy neutrinos goes through the star. In what fol- 
lows we characterize the general physical situation and es- 
timate the gamma-ray and neutrino output of these sys- 
tems, under a variety of conditions. The results will be 
of interest for current ground-based TeV gamma-ray tele- 
scopes and future satellites for MeV-GeV gamma-ray as- 
tronomy, as well as for km-scale neutrino detectors. 



2. Jet-wind interaction 

In order to make specific calculations we will consider a 
misaligned MQ with a high-mass companion. The jet axis 
is assumed, for simplicity, to be parallel to the orbital 
radius a, and pointing toward the massive star of radius 
R+. As in Romero et al. (2003), we consider an expanding 
conical jet, i.e. the radius of the cross-section is Rj(d) — 
dd/2, where 9 is the opening angle of the jet (in radians) , d 
is the distance from the compact object, and r is the radial 
coordinate with origin at the center of the star (d = a — r). 
Assuming adiabatic expansion, the magnetic field inside 
the conical jet changes with d as 



B- J {d)=B- i {d )[ -f 



(1) 



where do is the distance from the jet's apex to the com- 
pact object. We will assume d — 50i? g , being R g the 
gravitational radius of the accreting star. 

In order to determine the kinetic power of the jet we 
will adopt the jet-disk coupling hypothesis proposed by 
Falcke & Biermann (1995) and applied with success to 
AGNs, i.e. the total jet power scales with the accretion 
rate as Lj = gjAfdiskC 2 - We will assume that relativis- 
ts protons provide the dominant contribution to the jet 
power, i.e. L 3 w 7buikAp"ipC 2 , where N p is the proton in- 
jection rate, m p the proton mass, and Tbulk is the bulk 
Lorentz factor of the flow. The presence of relativistic 
hadrons has been inferred from iron X-ray line observa- 
tions for the case of SS433 (e.g. Kotani et al. 1994, 1996; 
Migliari et al. 2002). Although direct and clear evidence 
exists only for this source so far, several models assume 
a hadronic content for the jets of some MQs (see Romero 
2004 and references therein). 

The high-mass companion star loses a significant frac- 
tion of its mass through a very strong supersonic wind, and 
then the mass transfer occurs mainly through this stellar 
wind. To ensure that there is no Roche lobe overflow the 
parameters of the binary system must satisfy the condi- 
tion R L > R+, where the Roche radius can be expressed 
in good approximation as 



R L ^a [0.38 + 0.2 log(M*/M co )] 



(2) 



for 0.3 < (M*/M co ) < 20 (Kopal 1959). Here, M* and 
M co are the masses of the primary massive star and the 
compact object, respectively. 

Typical mass loss rates and terminal wind velocities 
for early-type stars are of the order of 10~ 5 M Q yr _1 and 



2500 km s , respectively (Lamers & Cassinelli 1999). The 
matter field will be determined by the continuity equation: 



M* = Anr 2 p(r)v(r), 



(3) 



where p is the density of the wind and v(r) = t>oo(l — 
R*/r)" is its velocity; i>oo is the terminal wind velocity, 
and the parameter f3 is ~ 1 for very massive stars (Lamers 
& Cassinelli 1999). Hence, assuming a gas dominated by 
protons, the particle density results: 



n(r) 



A.-Km v r 1 v(r) 



(4) 



The pressure exerted by the stellar wind at a distance 
r from the star is then: 



-fwind — 



(5) 



47rr 2 ' 

whereas the pressure due to the magnetic field of the star 
is 



Pmag = B(r) 2 /87T, 



(6) 



where B{r) = ^JB^ + B 2 . The magnetic field of the 

star can be represented by the standard magnetic ro- 
tor model (Weber and Davis 1967; White 1985; Lamers 
and Cassinelli 1999): B^/Br — (u*/woo)(l + r /R*) and 
B r = B+(R+/r) 2 , where v± = (0.1 — 0.2)^00 is the rota- 
tional velocity at the surface of the star (Conti & Ebbets 
1977), and B* is the corresponding stellar magnetic field. 

If the jet propagates with an opening angle 9, at a 
distance d from the compact object the ram pressure of 
the jet plasma is 



Pi 



u 



:9 2 d 2 ' 



(7) 



As shown by Bednarek & Protheroe (1997) for the case 
of AGN's jets, a strong stellar wind can stop the jet flow 
before it can reach the stellar surface. The interaction re- 
gion will be formed by a contact discontinuity and two 
quasi-parallel shocks, one propagating backwards through 
the jet and the other moving toward the star through the 
wind. 

The position r s where the shock front will form can be 
obtained by equating the expression for the jet pressure to 
the expressions given in Eqs. (0 or (JSJ). The location of the 
double-shock structure will be determined by the specific 
parameters assumed for the jet and the stellar wind. As an 
example, we will work out the details of a family of models 
sharing the parameters listed in Tabled For these values, 
it happens that Rl ~ 1.53 i?* and there is no Roche lobe 
overflow. 

With this set of basic parameters, the ram pressure of 
the jet can not be balanced by the wind when <j > 0.006. 
Above this value, however, the magnetic pressure can still 
balance and stop the jet. For gj = 10~ 2 , a magnetic field of 
B+ > 200 G is necessary to form a terminal shock front for 
the jet in the wind region. Early-type stars can support 
surface magnetic fields of even ~ 10 3 G as reported by 
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Donati et al. (2002). In order to study systems where the 
jet is balanced by the wind, we will assume different val- 
ues for the stellar magnetic field and the jet power. After 
discussing these systems we will turn to MQs with very 
strong jets in Section 

Table 1. Basic parameters assumed for the model 



Parameter 


Symbol 


Value 


Mass of the compact object 


M co 


10 M Q 


Mass of the donor star 


M* 


46.2 M Q 


Jet's injection point 


do 


50 Rl 


Opening angle of the jet 





5° 


Bulk Lorentz factor 


'ytv.ik 


10 


Radius of the companion star 


R* 


15 R @ 


Stellar effective temperature 


T e g 


40000 K 


Mass loss rate 


M* 


10" 5 M Q yr" 1 


Terminal wind velocity 




2500 km s _1 


Compact object accretion rate 


Mdisk 


10" 8 M Q yr" 1 


Wind velocity index 




1 


Efficiency of the shock 


X 


0.1 


Orbital radius 


a 


3 R* 


L R g = GM co /c z . 



3. Particle re-acceleration and diffusion 

The shocks generated by the interaction between the rel- 
ativistic plasma in the jet and the stellar wind can accel- 
erate particles by Fermi mechanism. For acceleration at a 
parallel shock we can write the acceleration rate as (e.g. 
Biermann & Strittmatter 1987, Protheroe 1998): 



dE 



= £ZecB, with £ ~ 0.015(u/c) 2 , 



(8) 



where u is the shock velocity and all values are in cgs 
units. 

The highest energy gained by the particles in the pro- 
cess is imposed by the different energy losses (synchrotron 
emission, pp interactions and photopion production) or by 
the available space in the acceleration region. For the pa- 
rameters considered in this paper the dominant energy 
loss is due to pp interactions with the stellar wind mat- 
ter that pervades the acceleration region. For a relativistic 
proton (E > 1.22 GeV) this energy loss can be estimated 
as in Mannheim & Schlickeiser (1994): 



dE 



0.65cn(r s )a pp (E — m p c 2 ), 



(9) 



where n(r s ) is the particle density of the wind close to the 
shock front, and a pp is the inelastic cross section for pp 
interaction at the energy E. The maximum energy for the 
protons then results from 



min[1.9xl0 15 ^y5#)eV 
\cJ G 

3xl0 14^o) eV] 



(10) 



We shall adopt u/c = 0.1 for the sub- relativistic shock 
and u/c — 0.5 for a relativistic case. For the former i?™ ax 
is determined by the losses, whereas in the second case it 
is fixed by the size constrain, i.e. the maximum gyroradius 
for a proton is ~ (a — i?*). 

The resulting proton spectrum from the diffusive ac- 
celeration will be a power-law extending from a minimum 
energy (E™ a - 10 GeV 1 ) up to E™ x : 



N p {E) = K T E 



-r 



10 m p c' 



<E< £™ ax . 



(11) 



The index r is ~ 2 for a strong non-relativistic shock and 
~ 1.5 for a relativistic shock (e.g. Protheroe 1998). 

Some fraction x of the energy transported by the 
jet will be transferred to the particles accelerated at the 
shock front. We will assume here an efficiency of x = 0.1 
(Blandford & Ostriker 1978). Then, the constant K T in 
(|llfl can be obtained from 



Xh = (a - R*fc I N P {E)E dE. 

'10m„c 2 



(12) 



The resulting population of relativistic protons will dif- 
fuse through the stellar wind, but not all of them will 
penetrate up to the densest parts where the probability 
of pp interactions is higher because of the modulation ef- 
fect of the wind (Romero & Torres 2003, Torres et al. 
2004). Whether the particles can reach the base of the 
wind or not is given by the ratio (e) of the diffusion and 
convection timescales, which are given by t d = 3r 2 /D and 
t c = 3r/v(r), respectively. Here D is the diffusion coeffi- 
cient defined as (Ginzburg & Syrovatskii 1964): 



D 



-Ai»c, 



(13) 



where A r is the mean-free-path for diffusion in the radial 
direction. Only particles for which e < 1 will be able to 
overcome convection and enter the dense wind region to 
produce 7-rays through pp interactions. 

The value of the diffusion coefficient is of course not 
well-known, especially for such a complex and turbulent 
medium as the wind of a massive star. Following White 
(1985) we will make use of the approximation given by 
Volk and Forman (1982), assuming that the mean- free- 
path for scattering parallel to the magnetic field direction 



is A 



10r e 



lOE/eB, where r g is the proton gyro- 



radius and E the proton energy. The value of A in the 
perpendicular direction is shorter: Aj^ ~ r g . The mean- 
free-path in the radial direction is then (sec Torres et al. 
2004): 



A r = r 9 (lOcos 2 6> + sin 2 0), 
where 

cos- 2 9= 1 + {B^/Brf. 



(14) 



(15) 



1 From Table 1 the Lorentz factor of the jet is ~ 10 and then 
those protons at the lower energy tail of the co-moving particle 
energy distribution (for which the microscopic Lorentz factor 
is ~ 1) are injected at the shock with energies ~ 10 GeV. 
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Through e = 1 we find now the minimum energy 
E™ m (r) below which the particles are taken away by con- 
vection in the wind (see the explicit expression for e in the 
paper by Torres et al. 2004). In what follows, we will take 
E p lm (r s ) as the minimum energy of a proton accelerated 
at the shock front to collide with a proton of the inner 
wind. 

4. Gamma-ray production 

We can now proceed to calculate the ir° gamma-ray lumi- 
nosity from the inner wind. The differential 7-ray emissiv- 
ity from 7r°-decays can be approximated by 

2Z {T) 

<Z 7 (£ 7 ) = A-ku w {E p )^^ J p (E 7 ) VA e(E p -E™ in ), (16) 

at the energies of interest here. Z^J^o is the so-called spec- 
trum weighted moment, defined as: 

Z JZ*° = / x r - 1 F p7T o(x)dx, (17) 
Jo 

where x = E x o / E p and F p7! o (x) is the dimensionless inclu- 
sive cross section (Gaisser 1990, p. 31) . The parameter tja 
reflects the contribution from different nuclei in the wind 
(for the standard composition of an early-type star we will 
assume t/a ~ 1.5). J p (E 1 ) = (c/Atv)N p (E 1 ) is the proton 
flux distribution evaluated at E — E 7 (units of protons 
per unit time, solid angle, energy-band, and area). 

The inelastic cross section for pp interactions at en- 
ergy E p sa Q^oE^/K, where K ~ 0.5 is the inelasticity 
coefficient and ^.0 = l.l(£ , p /GeV) 1 / 4 is the multiplicity 
for the neutral pions (e.g. Ginzburg & Syrovatskii, 1964), 
can be represented above E p w 10 GeV by the expres- 
sion given in Aharonian & Atoyan (1996). Typically, at 1 
TeV, a pp ~ 34 mb. Finally, Q(E p - £™ in ) is a Heaviside 
function that takes into account the fact that only pro- 
tons with energies higher than £?™ m (r s ) will diffuse into 
the densest part of the wind. The spectral 7-ray intensity 
(photons per unit time and energy-band) is: 

Iy(E y ) = J n{r)q^(E^)dV, (18) 

where V is the interaction volume. The generated lumi- 
nosity in a given band (E2, E±) then is 

rE 2 

L 7 = / E 1 I 1 (E 1 )dE 1 (19) 

(see, e.g., Romero et al. 2003 for details). 

In order to estimate the 7-ray luminosity that could 
be produced in the jet-wind interaction of a misaligned 
MQ with a high-mass companion, we have fixed the pa- 
rameters that describe the geometry of this scenario for 
a standard case (Tabic . Note that the radius of the 
star holds a ~ 26°. 5 angle as seen from the compact ob- 
ject, so the star fills ~ 10.5% of the compact object's sky. 
We have calculated the luminosity in high-energy photons 



for different jet-disk coupling values (gj), and for differ- 
ent values of B+ and Bj(do). The most effective particle 
acceleration occurs at the relativistic plasma. We have 
considered a sub-relativistic strong shock that leads to 
a spectrum with r ~ 2 and a mildly relativistic shock 
with u — 0.5c (r ~ 1.5). Table |21 characterizes (in the 
first three columns) a number of specific models, ordered 
by decreasing jet power, and then gives the location of 
the shock front, the maximum and minimum energy for 
protons penetrating into the base of the wind, and the 
integrated gamma-ray luminosities generated for the sub- 
relativistic and the relativistic cases of particle accelera- 
tion. 

The spectral high-energy distributions for some of the 
computed models are shown in Figure ^ The plots corre- 
spond to L 1 (E 1 ) — E 2 I 1 (E 1 ) 1 with an exponential- type 

cutoff at E y ~ (^ ax /GeV) 3 / 4 /20- The obtained isotropic 
luminosities above 1 TeV are in the range 10 30 — 10 32 erg 

High-energy gamma rays can be absorbed in the stel- 
lar photon field through pair production, and then inverse 
Compton emission from these pairs can initiate an e ± -pair 
cascade. This effect has been studied in detail by Bednarek 
(1997, 2000) and more recently by Sierpowska & Bednarek 
(2005). Figure [3 shows the r 7 opacity of the stellar pho- 
tosphere to gamma-ray propagation as a function of the 
gamma-photon energy and the distance r from the star 
where it is created. This opacity, for a photon moving ra- 
dially from the massive star, can be roughly computed, 
averaging over the injection angles, as: 

/>oo />oo 

T 1 (E 1 ,r) = J J N(E„ r')<j e - e+ (E*,E 7 )dE*dr', (20) 

where E+ is the energy of the photons emitted by the star, 
E.y is the energy of the 7-ray, and a e - E +(E ir ,E 1 ) is the 
cross section for photon-photon pair production (e.g. Lang 
1980). The stellar photon distribution is that of a black- 
body peaking at the star's effective temperature (T e g): 

N(E*, r) = {nB(E+)/hE*c)R 2 Jr 2 , (21) 

where h is the Planck constant and 

£?(£*) = [2£//(/ic) 2 ]/(e B * /fcToff - 1). (22) 

The inverse Compton opacity Tic(E e , r) determines 
the advance of the cascades after the pair creation. We 
have estimated that gamma-rays with energies below ~ 
10 12 eV will be mostly absorbed, triggering IC cascades. 
More energetic photons (> 10 TeV) will escape. The re- 
sult will be a soft spectrum at GeV energies (Bednarek & 
Protheroe 1997) and a peak at high energies (Figure 0). 

5. Neutrino production and possible detection 

The v^ + v^ neutrino intensity I y {E v ) may be derived from 
the spectral 7-ray intensity by imposing energy conserva- 
tion (Alvarez-Muhiz & Halzen 2002, Torres et al. 2004): 

J E^E^dE^ = C J ' E v I v (E u )dE Vl (23) 
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Table 2. Some results about the high-energy photon emission in misaligned MQs. The A, B, C letters naming the 
model indicate a high, middle or low jet-disk coupling, respectively. 



Model 


<h 








r s 


Trunin 


irim ax 
p, sub 


771 max 

"p, rel 


L -y,(T=2) 


^ 7 ,(r=i.5) 








[G] 


[G] 


[Re] 


[TeV] 


[TeV] 


[TeV] 


[ergs -1 ] 


[ergs -1 ] 


A 1 


1U 


12 


ouu 


i n 


91 7 


Q9 9 
yz . z 


i an 


OUUU 


991 n 32 

Z.Z 1U 


a i in 32 




i n~ 

1U 


2 


ouu 


i nn 

1UU 




Q9 9 
yz.z 


i onn 
iyuu 


OUUUU 


7^1 n 32 

( .0 1U 


q i n 32 


T> 1 


1U 


3 


i 

1 


i n 


Q A 9 


n o 
U.z 


ryu 


oUUU 


1.9 1U 


a i n 3 l 
0.4 1U 


132 


10" 


3 


1 


100 


34.2 


0.2 


1900 


30000 


2.1 10 32 


4.3 10 31 


B3 


10" 


3 


100 


10 


34.2 


21.3 


190 


3000 


6.4 10 31 


4.7 10 31 


B4 


10" 


3 


100 


100 


34.2 


21.3 


1900 


30000 


1.1 10 32 


4.1 10 31 


C'l 


10" 


4 


1 


10 


41.1 


0.2 


190 


3000 


1.9 10 31 


5.4 10 30 


C2 


10" 


4 


1 


100 


41.1 


0.2 


1900 


30000 


2.1 10 31 


4.3 10 30 


C3 


10" 


4 


100 


10 


41.1 


20.8 


190 


3000 


6.5 10 30 


4.7 10 30 


(74 


10" 


4 


100 


100 


41.1 


20.8 


1900 


30000 


1.1 10 31 


4.1 10 30 




8 10 12 14 8 10 12 14 

Log E [eV] Log E [eV] 



Fig. 1. Spectral high-energy distribution at TeV energies from the jet-wind interaction in a misaligned MQ. The left 
panel corresponds to T = 2, and the right to T = 1.5 (see text). The dashed lines correspond to the production spectra 
whereas the solid lines represent the absorbed spectra. 



where the limits of the integrals are (-E^K), the 

minimum (maximum) energy of the photons (neutrinos) 
and C is a numerical constant of order one (i.e. the to- 
tal neutrino luminosity produced by the decay of charged 
pions is a comparable quantity with that of the photons 
Lj). The neutrino intensity will be given by 

UE V ) = C v (E v /evg)- r e~ E »l E ™\ (24) 

where the constant C v is in erg -1 s _1 . The values of 
C„ (r=i.5) and C J/i (r=2) f° r the different models are listed 
in Table 

The neutrinos are produced in the inner wind around 
the star. A fraction of them, which will depend on the 
specific geometry of the system, will be intercepted by the 
star. The star, in turn, can absorb some of these neutri- 
nos. We have estimated the neutrino energy deposition 
rate r\ v (see Gaiser et al. 1986) inside the star. The results 



confirm that no effects of neutrino heating will be observ- 
able during the stellar lifetime. This fact is expected since 
only a small fraction of the power of the jet is directed 
toward the star by neutrinos {y\ v ~ 10 27 — 10 30 erg s^ 1 ) 
and this represents a very small fraction (10 -11 — 1CP 8 ) 
of the luminosity of the star. 

Although the neutrino emission has no effect upon the 
companion star, perhaps it could be detectable at the 
Earth. If the distance to the system is c^mq, the neutrino 
flux arriving to a km-scale detector like ICECUBE will be 
F V {E V ) — I v (E L ,)/(4nd^ l Q). To make some quantitative 
estimate we adopt in what follows g?mq ~ 5 kpc. 

The neutrino signal can be estimated as (Anchordoqui 
et al. 2003) 

S = T ohs [ dE v A ef! F v {E u )P v ^(E v ), (25) 
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r ~~~ - 




- 2 R, - 
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9 


10 11 


12 


13 14 



LogE y [ eV] 

Fig. 2. Stellar photosphere's opacity due to photo-pair creation, as function of the 7-ray energy and the distance from 
the star, for a photon moving radially from the massive star. Upper-right: Optical depth for fixed values of energy. 
Lower panel: Optical depth for fixed distances. 



where P^^E^) » 1.3 x 1CT 6 (E„/TeV) denotes the 
probability that a v of energy E v ~ 1 — 10 3 TeV, on a tra- 
jectory through the detector, produces a muon (Gaisser et 
al. 1995). T bs is the observing time and A e g the effective 
area of the detector. On the other hand, the noise will be 
given by 

1/2 



N 



T obs / dE v A en Fn{E v )P v ^{E v )m 



(26) 



where Afi is the solid angle of the search bin (Af2io xl o ss 
3 x 1CT 4 sr for an instrument like ICECUBE, Karle 2002) 
and 



F B (E V ) < 0.2(£„/GeV) 



-3.21 



GeV" 



cm 



sr 



(27) 



is the z/ M + atmospheric i/-flux (Volkova 1980, Lipari 
1993). Using as the integration limits the minimum and 
maximum energies detectable by the instrument (i.e. 1-10 5 
TeV) we obtain T Q b s for a signal-to-noise ratio S/N = 3. 
The results are in all the cases (models Al to C4) greater 
than several thousands years, rendering the detection im- 
practical with the current sensitivity. Once again, these 
are lower limits, since the duty cycle is not 1. Only in the 
case of very powerful jets with > 0.1 a detection might 
be possible on human-life timescales. We will discuss this 
case in the next section. 



Table 3. Some results about the neutrino emission in mis- 
aligned MQs. Case Dl corresponds to the direct impact 
of a jet with = 0.1 onto the stellar surface. 



Model C 



a 





erg x s 1 


erg 1 s 1 


Al 


4.0 10 32 


1.7 10 32 


A2 


3.3 10 32 


5.4 10 31 


Bl 


3.7 10 31 


1.6 10 31 


B2 


3.1 10 31 


5.4 10 30 


B3 


3.9 10 31 


1.7 10 31 


B\ 


3.3 10 31 


5.4 10 30 


CI 


3.7 10 30 


1.6 10 30 


C2 


3.1 10 30 


5.4 10 29 


C3 


3.9 10 30 


1.7 10 30 


C4 


3.3 10 30 


5.5 10 29 


Dl 


3.7 10 32 


3.1 10 34 
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6. Powerful jets and direct interaction with the 
star 

We consider now a heavy jet directly impacting the star 
(qj = 0.1 2 ). The proton spectrum of the relativistic jet flow 
is assumed to be a power law N' p (E' p ) = K T E' p , valid for 
mm < _£/ < E p max , in the co- moving jet frame. The cor- 
responding particle flux will be J' p (E' ) = (c/Att)N p (E' ). 
Since the jet expands (in a conical way), the proton flux 
can be written as: 



W - 5*" (t 



EL 



(28) 



where n is positive and we assume n — 2 (which corre- 
sponds to the conservation of the number of particles, see 
Ghiscllini et al. 1985), arid a prime refers to the jet frame. 
Using relativistic invariants, it can be proved that the pro- 
ton flux reaching the surface of the star, in the lab frame, 
becomes (e.g. Romero et al. 2003) 



Jp{,Ep 



cK r 

47T 



d 



-R* 

7bulk (E p + A, 



1 



/3b E p 



(29) 



where 7b u ik is the jet Lorentz factor, /3b is the correspond- 
ing velocity in units of c, and we have considered protons 
with small transverse momentum (hence the lack of de- 
pendency in the particle angle when compared with the 
corresponding expression in Romero et al. 2003). 

The number density uq' of particles flowing in the jet 
at Rq = R(do) can be determined as in Romero et al. 
(2003) and the normalization constant Kr can be derived 
provided that £'™ ax » E'™ in . Since E'f n ~ 1 GeV, this 
is always the case, and 



/mm\ p — 1 



K T = n '(T - l)(E>™ n ) 



(30) 



In the numerical calculations we have considered _B'™ ax = 
100 TeV, 7buik = 10, and two cases for the proton energy 
distribution in the jet 3 : T = 2 and T = 1.5. 

The gamma-rays generated by the jet-star interaction 
will be absorbed by the electromagnetic cascades in the 
star. The main signature of such cascades could be a rather 
broad e + — e~ annihilation line. The neutrino flux from 
charged pion decays will be significantly larger than in the 
case of the jet-wind interactions. For the density profile of 
the star we have used 



p(r) =p c [l- {r/R+f 



(31) 



2 Lj = q s 5.7 10 38 erg s" 1 

3 Notice that these are the original particle distributions in 
the jet and not the result of re-acceleration at the colliding 
shock region as before. Here the jet impacts directly in a high- 
density region. 



with p c ~ 0.05 g/cm 3 , in accordance with the adopted 
stellar mass. 

In the lab frame, the neutrino energies reach ~ 1 TeV, 
then only the hard energy tail of the neutrino spectrum 
lies in the range of observation of ICECUBE. The neutrino 
energy deposition rate for T — 1.5 results r/^ ~ 2.4 1 33 erg 
s _1 . The neutrino signal on the Earth produces a number 
of muons per year of = 3.25. This means a detectable 
signal at 3cr- level for an observing time T D b s ~ 15 years, 
in the case of a source at 2 kpc and a close alignment of 
the jet with the line of sight and a duty cycle of 20%. 

7. Comments 

The type of misaligned MQ discussed in Section 01 is a 
variable gamma-ray source. For the masses and orbital 
radius given in Tabled the orbital period is ~ 2.5 days. 
Since the star occupies ~ 10 % of the sky of the compact 
object and there are two jets in the orbital plane, the 
expected duty cycle is ~ 20 %. The gamma-ray emission 
is ignited twice per orbit, every ~ 30 hours. The source 
remains active for around ~ 6 hours each time, since the 
diffusion timescale for the accelerated particles is much 
shorter than this value: t d (S™ ax ) < 29 s and t d (£'™ in ) & 
203 s for the models considered here. All these numbers 
will change, of course, from one case to another depending 
on the specific configuration, but in all cases the result 
will be a variable source unless the diffusion timescale is 
much longer than the interval between the different jet- 
wind collisions. 

Because of the isotropization of the particle flux at the 
shock front located at r s the gamma-ray emission will not 
be as much dependent on the viewing angle as in the case 
of pure relativistic jet models (e.g. Romero et al 2003). 
Occultation events will occur only if the observer lies close 
to the orbital plane. Such events would be rare since they 
require an alignment among compact object, star and ob- 
server. The flux, additionally, would not be completely 
suppressed since the relativistic particles diffuse around 
the entire inner wind. 

Since misaligned MQs can be TeV and (through elec- 
tromagnetic cascades in the stellar photosphere) MeV- 
GeV sources, we could ask whether some of the unidenti- 
fied gamma-ray sources detected by EGRET instrument 
might be associated with them. There are ~ 20 unidenti- 
fied variable sources at low galactic latitudes in the Third 
EGRET Catalog (Hartman et al. 1999, Torres et al. 2001, 
Nolan et al. 2003). Bosch- Ramon et al. (2004) have shown 
that they form a distinct population with a distribution 
consistent with what is expected for high-mass MQs. The 
total number of high-mass X-ray binaries brighter than 
10 34 erg s _1 at X-rays in the Galaxy is estimated in ~ 380 
(Grimm et al. 2002). Only a small fraction of them will 
have the type of persistent jets necessary to produce de- 
tectable gamma-ray sources. Around 20 sources seems to 
be a reasonable estimate (Bosch- Ramon et al. 2004). If 
the original orientation angles of the jets are randomly 
distributed we could expect that in 10-20 % of these sys- 
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terns the jet could be periodically directed towards the 
companion star (Butt et al. 2003). This means that 2-4 
of the variable EGRET sources at low galactic latitudes 
might be associated with misaligned MQs. Hints of peri- 
odic variability could help to identify these sources. Their 
main footprints, however, would be a stronger TeV emis- 
sion than what is expected in normal MQs. Hence, fu- 
ture observations with Cherenkov telescopes like HESS, 
MAGIC or Veritas might lead to the identification of these 
sources. 

8. Conclusions 

The jets of MQs can form a small angle respect to the or- 
bital plane. In some high-mass X-ray binary systems the 
jet could be directed straight to the companion star. We 
have shown that in such a situation the wind of the star 
can balance the jet pressure in many cases. A jet-wind 
collision region is formed between the star and the com- 
pact object, with a contact discontinuity separating both 
shocked media. Relativistic particles injected by the jet 
in the strong shock can be re-accelerated up to very high 
energies and isotropized in this region. These particles dif- 
fuse into the inner wind as far as diffusion dominates over 
convection. Gamma-rays and neutrinos can be then pro- 
duced in the densest part of the wind. We have presented 
calculations of the gamma-ray luminosity, which might 
be detectable through high-energy telescopes. The neu- 
trino flux is usually below the sensitivity of km-scale de- 
tectors like ICECUBE unless extremely powerful jets are 
assumed. Such jets could reach the star, triggering cas- 
cades and producing nucleosynthesis there, as discussed 
by Butt et al. (2003). We have also shown that photon- 
photon absorption of gamma-rays with energies below 10 
TeV can initiate electromagnetic cascades in the stellar 
photosphere, leading to MeV-GeV sources that might con- 
tribute to the population of variable gamma-ray sources 
detected by EGRET at low galactic latitudes. Future in- 
struments like GLAST will help to unveil the nature of 
these sources. 
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